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Abstract 

We study various time series of surface layer wind velocity at different locations and provide 
evidences for the intermittent nature of the wind fluctuations in mesoscale range. By means of the 
magnitude covariance analysis, which is shown to be a more efficient tool to study intermittency 
than classical scaling analysis, we find that all wind series exhibit similar features than those 
observed for laboratory turbulence. Our findings suggest the existence of a "universal" cascade 
mechanism associated with the energy transfer between synoptic motions and turbulent microscales 
in the atmospheric boundary layer. 
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I. INTRODUCTION 



Atmospheric surface layer motions are a source of many challenging problems. The issue 
of designing a faithful statistical model of spatio-temporal wind speed fluctuations has been 
addressed in various fields like boundary layer turbulence phenomenology, meteorology, wind 
power control and prediction or climatology. From turbulent gusts to hurricanes, breezes 
to geostrophic wind, the wind process is characterized by a wide range of spatio-temporal 
scales and all the above mentioned disciplines mainly focus on a specific subrange of scales. 
The modeling of wind speed behavior in the mesoscale range is of great interest for example 
in wind power generation or in order to control polluant dispersion. In this range of scales 
extending from few minutes to few days (~ 1-1000 km), the properties of wind velocities are 
less known than in the range of planetary motions (synoptic scales) or turbulent motions 
(microscales) . From a physical point of view, because of the importance of boundary 
conditions, the heterogeneous and non-stationary nature of the processes involved, it is well 
admitted that mesoscale wind regimes strongly depend upon various factors like atmospheric 
conditions, the nature of the terrain and may involve periodic variations (caused by diurnal 
temperature variations). Unlike microscale Kolmogorov homogeneous turbulence, mesoscale 
fluctuations are therefore not expected to possess any degree of universality However, 
during past few years, some papers have been devoted to the analysis of scaling laws and 
intermittency features at large scales in many geophysical signals like temperature, rainfall 
or wind speeds 0, O, Q|. In Ref. Q], the authors showed that surface layer wind velocities 
recorded at low frequency using a cup anemometer display multiscaling properties very much 
like in the high frequency turbulent regime. Moreover, they claimed that random cascade 
models could be pertinent to reproduce the observed intermittent fluctuations. Along the 
same line, in a Multifractal Detrended Analysis was performed on four different hourly 
wind speed records and revealed some multiscaling properties of the series. Even if these 
studies did not go beyond simple scaling exponent estimation and did not consider problems 
related to the statistical significance of the obtained results, they had the merit to address 
questions about the possible intermittent nature of wind variations in the mesoscale range. 
As reviewed below, one of the consequences of multiscaling and intermittency is that small 
scales fluctuations are strongly non Gaussian and characterized by "bursty" behavior. In 
Refs. jsl, 9], such features have been precisely observed on wind variations statistics at largest 
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microscales (or at smallest mesoscales) and have been shown to be related to a "fluctuating" 
log-normal turbulent intensity. 

In this paper we suggest that this "fluctuating turbulent intensity" results from a random 
cascading process initiated at a time scale of few days. Our aim is to show that, in some 
sense, "turbulent" cascade models are likely to be pertinent at larger scales, in the so-called 
mesoscale regime. As compared to the previously cited papers, our analysis relies upon the 
use of magnitude (i.e. logarithms of velocity increments amplitudes) correlation functions. 
From a mathematical point of view, long-range correlated magnitudes have been shown to be 



at the heart of the construction of continuous cascade processes 
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ll| . For a practical point 



of view, magnitude covariance possesses interesting properties as far as statistical estimation 



problems are concerned 
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14| . We show that these correlation functions can be reliably 



estimated and are very similar to those associated with longitudinal velocity time series of 
laboratory turbulent experiments. Moreover we observe some universal features among the 
various analyzed series. 

The paper is organized as follows: in section [Til we make a brief review of intermittency 
and the related notion of random cascade models. We emphasize on the interest of studying 
magnitude correlations and discuss its relationship with multiscaling properties. We then 
compare, on an empirical ground, the relative performances of intermittency estimators 
relying upon scaling and magnitude covariance. This section ends with a review of a recent 
work of B. Castaing that has shown how intermittency in Lagrangian and Eulerian frames is 
"observed" on time series recorded at a fixed spatial position. Our main experimental results 
are then presented in section IIHI After a rapid description of various wind data we have 
studied, we show that wind surface layer variations in mesoscale range have intermittent 
properties and possess "universal" magnitude covariance similar to laboratory turbulent 
fluctuations. Discussions and prospects are reported in section [TV] 
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II. INTERMITTENCY: FROM MULTISCALING TO MAGNITUDE CORRELA- 
TIONS 



A. Multiscaling and intermittency 

Small scale intermittency is one of the most challenging problems in contemporary tur- 
bulence research. It is generally associated with two distinctive features: the first one is that 
at small scales the probability density functions (pdf) of velocity variations are strongly lep- 
tokurtic (with "stretched exponential" tails) while they are almost gaussian at larger scales. 
The second one is that the so-called structure functions display multiscaling properties. As 
we shall see, these two properties are in fact equivalent within the multiplicative cascade 
picture. Let us make a brief overview of these notions. 

We denote X(t) a continuous process (for instance the time variations of a velocity field 
component at a fixed location) and let S T X(t) be its increments over a scale r: 5 T X(t) = 
X(t + t) — X{t). One usually defines the order q structure function of X as 

S q {r) = J \8 T X(u)\idu (1) 

and the ((q) spectrum as the scaling exponent of S q (r): 

S q (r) ~ r«*> (2) 

If the function ((q) is non-linear one says that X(t) is a multifractal process or an intermit- 



tent process. In that shown e.g. in [15], £(g) is necessarily a concave function and 



the previous scaling holds in the range of small r ; r — > means precisely r«T where T is 
a coarse scale called the integral scale in turbulence (in general associated with the injection 
scale). The intermittency coefficient is a positive number that quantifies the non-linearity 
of ((q) and can be defined 1 as, e.g., 

A 2 = -C"(0) (3) 

The most common example of non linear ((q) function is the so-called log-normal spectrum 
which is a simple parabola: 

A 2 

C(?) =aq- —q 2 (4) 



1 Notice that one can find different definitions of the "intermittency coefficient" or the "intermittency 



exponent" in the literature (see 16| for example). However within the log- normal cascade model discussed 
below they are all equivalent. 
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In that case A 2 corresponds to the constant curvature of ((q). 

In order to estimate the multiscaling properties and/or the intermittency coefficient, 
one can directly estimate partition functions from the data and obtain the exponent ((q) 
from a least square fit of S q (r) in log- log representation. However this method suffers 
from various drawbacks. First, from a fundamental point of view, one has to distinguish 
temporal (or spatial) averages from ensemble averages. For instance, in the case of a log- 
normal multifractal, rigorously speaking, only the "ensemble" average (\S T X(u)\ q ) behaves 
as a power law with an exponent ((q) as given by Eq. (jlj). The temporal or spatial average 
has an exponent spectrum ((q) that becomes linear above some value of q and is no longer 



parabolic [17j. In order to estimate from moment scaling over a wide range of q one 
has to use the so-called "mixed" asymptotic framework [3] . But the main problem remains 
that high order moment estimates require very large sample size. Moreover the scaling 
behavior can also be altered by finite size effects, discreteness and non stationarity effects 
like periodic perturbations or periodic modulations of the data (see below). A more reliable 
method first introduced in |l2j (see also 19]) relies upon the so-called magnitude cumulant 
analysis. It simply consists in focusing on the scaling behavior of partition functions around 
q = 0. According to this approach, the structure function is written as: 

too £, \ 

where Ck{r) is the k-th cumulant associated with the random variable uj t (u) = 
ln(\S T X(u)\). The logarithm of the increment will henceforth be referred to as the mag- 
nitude of velocity increments. Note that C\ is simply the mean value of uj while Ci is its 
variance. Thanks to the scaling relationship ([2]), one deduces that all cumulants behave as: 

C k {j) = c fc ln(r) +r k (6) 

where the constants {r k } account for both the integral scale and the prefactors in the scaling 
relationship (TSJ). The function ((q) can therefore be expressed in terms of a cumulant 
expansion: 

k 

In particular one sees that the intermittency coefficient is directly involved in the behavior 
of the magnitude variance as: 

C 2 (r) = -A 2 ln(r)+r 2 . (8) 



Eq. (jHJ) has been successfully used to estimate the intermittency coefficient of longitudinal 
velocity fields in 3D fully developed regime under various experimental conditions. A com- 
mon value close to A 2 ~ 0.025 has been obtained. It is remarkable that this intermittency 



value appears to be universal 



12 



Let us end this brief review by discussing the relationship between intermittency and the 
small scale "bursty" behavior of velocity increments. Indeed, it is well known in turbulence 
laboratory experiments that large scale increment pdf are close to be normal while small 
scales have larger tails. The flateness strongly increases when one goes from large to small 
scales. We focus only on the log-normal CcLSG, 1.6., cL parabolic ({q), extension of our consid- 
erations to other laws being straightforward. Thanks to the structure function multiscaling 
interpreted as a moment equality, by simply performing a time scale contraction, r' = st 
(s < 1), one can write: 

5 ST X = e n °5 T X (9) 

Law 

where Q s is a Gaussian random variable of variance —A 2 ln(s) which law is denoted as 
G S (Q). By simply assuming that for r = T (T being the integral scale), the increments of X 
are nor mally distributed, one obtains the small scale pdf of 5 T X by the well known Castaing 
formula 21] : 



p(z,r) = (2tt)-* / G T/T (tt)e n e — dti (10) 

This means that, at scale r, the pdf of the increments of the process X(t) are obtained 
as a superposition of Gaussian distributions of variance e 2Q where Q is itself a gaussian 
random variable of variance increasing at small scales like A 2 ln(T/r). The smaller the scale, 
the larger the variance of Q and therefore the larger the tails of the increment pdf p(z,r). 
The continuous deformation from Gaussian at large scales towards stretched exponential 
like shapes at small scales, observed in laboratory turbulence experiments, has been shown 
to be well accounted by the transformation ffTO]) . 



B. Random cascades and logarithmic magnitude covariance 

A natural question that arises after the previous analysis is how can we explicitely build 
models that possess multiscaling properties? In other words, since multiscaling is equivalent 
to intermittency, how can the variance of the magnitude increase as a logarithm of the scale 
as described by Eq. (jSj) ? The answer comes from the self-similarity Eq. (Q that can be 
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iterated and interpreted as a random cascade: when one goes from large to small scales, one 
multiplies the process by a random variable W s = exp(Q s ). 

Usually one starts by building a non-decreasing (i.e. with positive variations) cascade 
process, denoted hereafter as M(t), which is often referred to as a multifractal measure 
though its variations are not bounded. More general multifractal processes (or multifractal 
"walks") can be simply obtained by considering a simple Brownian motion B(t) (or any self- 
similar random process) compounded with the measure M(t) considered as a stochastic time: 
X(t) = B [M(t)\. In finance literature M(t) is often referred to as the 'trading time' while 
in turbulence M(t) can be associated with the spatial or temporal distribution of energy 



dissipation. T 
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[see e.g. 
as models o 
Mandelbrot 
cascades 



re statistical properties of X(t) can be directly deduced from those of M(t) 



23|). Random multiplicative cascades measures were originally introduced 
the energy cascade in fully developed turbulence. After the early works of 



24, 
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26 L a lot of mathematical studies have been devoted to discrete random 



30, 



3l|. Let us summarize the main properties of these constructions, 



set some notations and see how these notions extend to continuous cascades. 

The simplest multifractal discrete cascades are the dyadic cascades defined by the follow- 
ing iterative rule: one starts with some interval of constant density and splits this interval in 
two equal parts. The density of the two sub-intervals is obtained by multiplying the original 
density by two (positive) independent random factors W of same law. This operation is then 
repeated ad infinitum. The integral scale corresponds to the size of the starting interval. A 
log- normal cascade corresponds to W = exp(f2) where Q is normally distributed. Peyriere 
and Kahane 27j] proved that this construction converges almost surely towards a stochastic 
non decreasing process M M provided (VTlnVT) < 1. The multifractality of (hereafter 
simply denoted as M) and of X(t) = B(M(t)) (B(t) being a Brownian motion) has been 
studied by many authors (see e.g. [29]) and it is straightforward to show that the spectrum 
of X(t) is 

C(«) = flf/2 - ln 2 {W^ 2 ) (11) 

The self-similarity Eq. ()§]) for dilation factors s = 2 n can also be directly deduced from the 
construction. 

Because the discrete cascade construction involves dyadic intervals and a 'top-bottom' 
construction it is far from being stationary. In order to get rid of theses drawbacks, con- 
tinuous cascade constructions have been recently proposed and studied on a mathematical 



ground 
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23 



32 



33j, |34|. They share exact multifractal scaling with discrete cascades 
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33j. With- 



but they display continuous scaling and possess stationary increments 
out entering into details, we just want to stress that these constructions rely upon a family 
of infinitely divisible random processes flj(t). The process M(t) is obtained as the weak 
limit, when I — > 0, of J*e Ql ^dv. In the log-normal case, fij is simply a Gaussian process 
defined by a covariance function that mimics the behavior of discrete cascade. Indeed, as 



shown in 



id, 



351 ] . if X(t) = B(M(t)) where M(t) is a discrete cascade as defined previ- 
ously with intermittency coefficient A 2 , then correlation function p(At) of the magnitude 
oj T {u) = ln(|<5 T A \u)\) decreases slowly as a logarithm function: 



p{At) d = Cov [lu t (u),lu t (u + At)] ~ -A 2 ln(— — ) for lags At < T. 



At 



(12) 



■T + t' 

The integral scale T where cascading process "starts" can therefore be interpreted as a 
correlation length for the variation log- amplitudes of X(t). The so-called Multifractal Ran- 
dom Walk (MRW) process introduced in Ref. [36] consists in constructing a log-normal 
multifractal process X(t) = B(M(t)) where B(t) is a Brownian motion and 



M(t) 



lim 

1^0+ 



3 Ci(M) du 



(13) 



where the fli(u) is a gaussian process with a logarithmic covariance as given in Eq. (fT2l . In 
Fig. Ufa) and[T](b) are plotted respectively a sample of a log-normal measure M(dt) (dt = 1) 
and a path of a MRW process X(t) corresponding to A 2 = 0.025 and T = 512. 

The studies devoted to continuous versions of discrete cascades have mainly shown that 
multifractal processes are related to exponential of logarithmic correlated random processes. 
The magnitude covariance function has been proven to be at the heart of the notion of 
"continuous cascade" . As we shall see below, it also allows one to estimate the intermittency 
coefficient in a more reliable way than methods based on scaling properties. 



C. Intermittency coefficient estimation issues 

As far as the problem of the estimation of the intermittency coefficient is concerned, it 
results from previous discussion that this coefficient, originally defined as the curvature of 
C(q) (Eq. (jHJ)), can be estimated either from the behavior of magnitude variance across 
scales (Eq. ([81) or from the slope of the time magnitude covariance in lin-log coordinates 
(Eq. m)- 
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FIG. 1: Estimation of the intermittency coefficient for a log- normal MRW multifractal measure (a) 
and a log-normal MRW process (b). In both cases the sample size is 32768 points, A 2 = 0.025 and 
T = 512. In (c) are plotted the magnitude covariance (•) and magnitude variance (o) as a function 
respectively of the logarithm of the lag and the logarithm of the scale. The slope of both curves 
provides an estimation of A 2 . One can see that for measure M(dt), the errors in the estimation 
are comparable. In (d) are displayed the same plots as in (c) but for the magnitude of the MRW 
process. Because of the additive noise and the small number of independent points at large scales, 
the estimation relying upon the magnitude variance turns out to be much more altered. 

Let us show that this latter definition is much more reliable than former one. The precise 
computation of the properties of estimators relying upon Eqs (JSJ) and (|12[) is a difficult 
task. In Ref. [371, it has been shown that a Generalized Method of Moments relying 
on the magnitude correlation function provides an unbiased and consistent estimator of 
A 2 . Estimators relying on Eq. (8]) have been precisely discussed within the context of 
atmospheric turbulence in Ref. 19J where the authors showed that it allowed one, by 
means of a bootstrap method called "surrogate analysis", to estimate A 2 and distinguish 
intermittent from non intermittent time series. 



A simple argument allows us to understand why magnitude covariance based estimation 
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(Eq. (IT21 ) is better than scale magnitude cumulant analysis (Eq. (jSJ)) in the case of mul- 
tifractal random walks. This is illustrated in Fig. [I] where both estimators are compared 
for a sample of a multifractal random measure constructed from a log-normal continuous 
cascade (Fig. [1(a)) and a MRW process (Fig. [1(b)). One sees in Fig. [He), that in the case 
one handles directly the random cascade measure, both estimators are roughly equivalent. 
However, in the case of a MRW process, the measure is not directly observed since it repre- 
sents the stochastic variance of the random walk. In this situation, the estimates based on 
magnitude covariance are better than those obtained from the magnitude variance scaling 
(Fig. [1(d)). This feature can be explained as follows: if one considers that for continu- 
ous cascades processes, Eq. flU]) is an equality in law for all finite dimensional distributions 
(f.d.d.), by taking the logarithm, one gets: 



where Q T (u) is a log-correlated gaussian random variable (Eq. (1121) ) which variance behaves 
like in Eq. (jSJ) while e(u) is a standardized normal white noise independent of Q T . Note 
that in the case of a multifractal measure, there is no noise term ln(|e|)) and the estimation 
error is small as compared to the situation of MRW when this term is present. One can 
then assume that the estimation error mainly comes from the presence of ln(|e|). This 
error is thus expected to decrease, in both CclSGS, clS -pFj= if N e is the number of effective 



independent samples of the considered signal. This effective number is proportional to the 
overall sample size but decreases with the scale r. Therefore, since at scale r, the effective 
number of independent points is N e ~ Nt" 1 , if one fits Eq. (JSJ) over the range of scale 
extending from r = 1 to r = r max , the mean error is roughly proportional to V'WxjV™ 1 - 
On the other hand, by fitting Eq. (Tl2|) at scale r = 1 from lags At = 1 to At = r max , one 
expects an approximately constant error of V A r_1 . Moreover, the error in the computation 
of the variance is proportional to the square root fourth moment of ln(e) while the error 
in the covariance is proportional to the second moment. Since the ratio between these two 
quantities is 3, the ratio of errors between the two methods is of order: 



u T (u) = \n(\S T X(u)\) 



f.d.d. 



n T (u)+\n(\e(u)\) 



(14) 



max 



(15) 



Let us notice that, rigorously, the error also depends on the overall range of values used 
to apply the linear regression and thus one expects, in the case of an estimation relying 
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upon Eq. (jSJ), there is an optimal range of scale to perform the fit. This optimum results, 
on the one hand, from a reduction of the error when one increases the fitting range of scales 
and, on the other hand, from the increasing of the error associated with the small number 



of independent points at large scales (see Ref. |38| where such a problem is considered in 
details). We have performed extensive Monte-Carlo estimations using synthetic log- normal 
continuous cascades and we have checked that for parameters, sample sizes and range of 
scales close to the experimental values considered below (see section ITO]) . the observed ratio 
between the errors is r exp ~ 10 in agreement with previous considerations. The intermit- 
tency coefficient of a random cascade is thus estimated in a much more reliable way using 
magnitude correlation functions. 



D. Squared log magnitude covariance as the result of Lagrangian and Eulerian 
intermittency 

In most fluid mechanics experiments, one usually records one or several components of 
the velocity field using an anemometer, so one has only access to the values of the field at 
a fixed spatial position as a function of time. In order to make inferences on the spatial 
properties of the velocity, i.e., the Eulerian field, one generally invokes the Taylor Frozen 
hypothesis or, when the turbulence rate is large, the Tenekes sweeping argument according 
to which small scale velocity fluctuations are mainly caused by Eulerian variations advected 
by large scale random motions [39j. However, as explained below, for some functions, single 
point measurements cannot be linked so easily to their Eulerian counterpart. Hereafter, 



we reproduce the argument of B. Castaing 



40 



ill 



in order to understand the shape of 



the single point magnitude covariance if one supposes that both Eulerian and Lagrangian 
velocity fluctuations are given by a continuous cascade as described previously. 

Let us denote Q(x,t) the magnitude at time t and position x. It is important to notice 
that Q(x,t) is a local field that does not depend on any spatial or time scale (for example 
in turbulence, Q(x,t) can be considered as the logarithm of the dissipation field or of the 
velocity increments at a scale smaller than the Kolmogorov dissipation scale). 

Our goal is to compute 

Cav[Q(x,t),n(x,t + At)] (16) 
as a function of the time lag At for a fixed value of x. 
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If one supposes that both Eulerian and Lagrangian fields are log-correlated, i.e., well 
described by a continuous cascade model, then 

Cav[n(x,t),Q(x + r,t)] = A 2 ln(— ^— - ) 

Cov[Q(x(t),t),n(x(t + At),t + At)} = p 2 ki(-^- — ) 



where At and r are time and space lags, r] e and r/i are small scale spatial and temporal cut- 
offs (i.e. the Kolmogorov scales in turbulence), L and T are spatial and temporal integral 
scales and A 2 , p 2 are Eulerian and Lagrangian intermittency coefficients. 

Let us suppose that the fluid particle at position x at time t + At was at position x' at 
time t. Thanks to the above covariance formula, one can write: 

n(x,t + At) = p!n(x\t) + e (17) 

where e is a random variable independent of Q(x',t) and: 



p 2 \n(- 



T 



» = -Vaff <18) 
But (at least from a statistical point of view) r — \x' — x\ — V At where V is a 'typical' 
velocity (the mean or r.m.s. velocity) so that, if one assumes that L ~ VT: 

Q(x',t) = p 2 n(x,t) +e' (19) 

where e' is a random variable independent of Q(x, t) and: 

P2 = Var(fi) (20) 
Then, if e and e' are uncorrelated, since the correlation coefficient between Q(x, t) and 

Q(x,t + At) is the product P1P2, we have 

Cov [Q(x, t), Q(x, t + At)] = Var(fi)pip 2 , (21) 

and by taking into account the fact that Var(Vt s ) = p 2 ln(T/^) = A 2 \n(L/r] e ) one obtains 
finally, for lags such that V At 3> r\ e and At 3> rjf. 

Cov [Q(x, t), t + At)] ~ ^ In 2 (^j = ^ In 2 (^j . (22) 
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If one plots vCov as a function of ln(r) one expects a straight line of slope: 

(23) 




In one knows the "Reynolds numbers" — and — , the intermittency coefficients can be 



estimated from the slope r as: 

A 2 = r 2 ln(^) (24) 
^ = r 2 \n(-) (25) 

Ve 

We see that, by taking into account both Eulerian and Lagrangian fluctuations on single 
point measurements, one should observe a squared logarithm magnitude covariance instead 
of the logarithmic behavior of Eq. ( 1T21) . This peculiar shape of magnitude covariance has 



indeed been precisely observed on laboratory fully developed turbulence data in Ref. 



12|. 



III. RESULTS: INTERMITTENCY IN MESOSCALE WIND FLUCTUATIONS 
A. The data series 

The results reported in the following are based on different wind velocity time series. The 
first data set consists in horizontal wind speeds and directions that have been recorded every 
minutes during 5 years (1998-2002) at our laboratory in Ajaccio - Vignola at a height of 10 
m by means of a cup anemometer. We also study hourly wind speed and direction data (10 
minutes averages) for 7 different sites in Corsica (France). The length of these series is also 
5 years. These data have been measured and collected by the French Meteorological Service 
of Climatology (Meteo-France) using a cup anemometer and wind vane at 10 meters above 
ground level. Finally, we also consider potential winds from KNMI HYDRA PROJECT 



available online [42j: they represent series of hourly (1 hour average), 10 meters potential 
wind speed gathered during several years at various locations in Netherlands. More specifi- 
cally we consider the series recorded at 3 different sites over a long period of several decades. 
Table [J summarizes the main characteristics of the studied series. 

In the sequel, v(t) will denote the modulus of the velocity horizontal vector while v x (t) 
and v y (t) will stand for its two components along arbitrary orthogonal axes x and y. We 
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Location 


Latitude 


Longitude 


Dates 


Sampling freq. 


Site 


Vignola (Ajaccio) 


41°56'N 


8°54'E 


1998-2003 


1 min 


50m, coastal, high hills 


Ajaccio 


41°55'N 


8°47'E 


2002-2006 


1 hour 


5m, coastal, plain, airport 


Bastia 


42°33'N 


9°29'E 


2002-2006 


1 hour 


10m, coastal, plain, airport 


Calvi 


42°31'N 


8°47'E 


2002-2006 


1 hour 


57m, coastal, hills 


Conca 


41°44'N 


9°20'E 


2002-2006 


1 hour 


225m, high hills 


Figari 


41°30'N 


9°06'E 


2002-2006 


1 hour 


21m, plain, airport, hills 


Renno 


42°11'N 


8°48'E 


2002-2006 


1 hour 


755m, mountains 


Sampolo 


41°56'N 


9°07'E 


2002-2006 


1 hour 


850m, mountains 


Eindhoven 


51°44'N 


5°41'E 


1960-1999 


1 hour 


20m, plain 


Ijmuiden 


52°46'N 


4°55'E 


1956-2001 


1 hour 


4m, coastal, plain 


Schipol 


52°33'N 


4°74'E 


1951-2001 


1 hour 


-4m, plain, airport 



TABLE I: Main features of the time series 




-5 -5 

ln(f) (hour -1 ) ln(f) (hour -1 ) 



FIG. 2: Power spectrum density of v x wind components 10 min rate Vignola series (solid line) and 
hourly Eindhoven series (dotted line) in log-log representation, (a) One clearly identifies the peaks 
associated with diurnal oscillation superimposed to an overall scaling regime where P(f) ~ 
with (3 ~ 1.6. This regime roughly extends from few days to few minutes. Turbulent atmospheric 
scales are below the finest resolved time scale, (b) Plot of the same spectra where the estimated 
local seasonal components have been removed. 

have by definition: 

v(t) = yjv x (t)* + v y (ty 

Because there is no well denned constant mean velocity direction with a small turbulent rate, 
we have chosen to study v x and v y separately instead of considering meaningless longitudinal 
and transverse components. 

The power spectrum analysis is one of the most common tools for analyzing random 
functions and is at the heart of a wide number of studies of wind velocity statistics. Since 
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V x (t) - 




v x (t) 
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FIG. 3: Time fluctuations of the v x {t) component of the Vignola velocity field, (a) At large time 
scales, one does not see any structure, we are in the flat (white noise) regime of the power spectrum, 
(b) At a finer time resolution, the series appears as a superposition of turbulent gusts and a more 
regular component. 



the pioneering work of Yen Der Hoven 



43|, 



44j , the shape of a typical atmospheric wind speed 



spectrum in the atmospheric boundary layer is still matter of debate. It is relatively well 
admitted that it possesses two regimes separated by low energy valley called the "spectral 
gap" located at frequencies around few minutes. This gap separates the microscale regime, 
where turbulent motions take place, from the mesoscale range. In the homogeneous turbulent 
regime, it is well known that the spectrum associated with the velocity field behaves like 
E(k) ~ fc -5 / 3 as predicted by Kolmogorov in 1941 [l| (k is the spatial wave-number). In 
the mesoscale range, for time scales greater than few minutes, the shape of this spectrum 
appears to depend on various factors. If some experiments suggest that a k~ 5 ^ 3 spectrum 



can extend up to synoptic scales 



surface layer 



45l . |46| in the free atmosphere, things are different in the 



471 ] . Some authors suggest that statistics a priori depend on local conditions 



(orographic, atmospheric,...) and one does not expect the same degree of universality as in 
the microscale regime. For example, in ref [3j, it is shown that the spectrum exponent may 
depend on the atmospheric stability conditions and also on the topography. 

In Fig. [2] are plotted the power spectrum of the v x component of two series: the "high 
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frequency" (10 min rate) series of Vignola and the hourly series recorded at Eindhoven. 
These two spectra that do not cover the same frequency range have been shifted by an 
arbitrary multiplicative factor. Series associated with the v y component or corresponding 
to other sites have similar features. One can see that, up to the mains peaks associated 
with diurnal wind oscillations (see below), these spectra are well described by a power-law 
P(f) ~ in a frequency domain which corresponds to time scales from few minutes to a 
characteristic time of few days (note that the low frequency behavior is much more reliable in 
the Eindhoven series since it covers a period close to ten time longer than the Vignola series). 
The value of the exponent is (3 ~ 1.6 for both sites. Within the framework of self-similar 
Gaussian processes |48|, the fact that (3 > 1 for all series indicates that the signal process 
v x (t) has continuous pathes. Since at low frequency the spectrum becomes flat, this means 
that the process appears regular at small scales and more "noisy" at larger scales. This 
feature is illustrated in Fig. [3] where we have plotted two samples of the component v x (t) of 
Vignola wind series over two different time intervals. At very large scale, the signal looks like 
a highly irregular "white noise" that corresponds to the power-spectrum low frequency flat 
behavior (Fig. H2(a)) while a zoom over a finer time interval reveals more regular variations 
(Fig. E(b)). Notice that one can also observe a daily oscillating behavior and high frequency 
turbulent gusts superimposed to these regular random variations. 

Since our goal is to study the stationary random components of the velocity field, we 
have preprocessed all time series in order to remove the seasonal components. In fact the 
diurnal oscillations S x>y are not really periodic but vary during the year according to the 
sun position. We have used a local harmonic parametrization of these components and 
computed the parameters by minimizing an exponential moving average of a quadratic error 
sec for more details). In order to study the fluctuations in these deseasonalized series 

and an eventual intermittency, we can compute, as in turbulence literature, its increments. 
Notice that one can alternaltively study velocity component wavelet coefficients or, if one 
wants to account for the low-frequency behavior, the error in the one step forward prediction 



of a Langevin like modeling of v x (t) and v y (t) [50(: 



V X ,y(t + 1) = S x>y (t + 1) + (1 - 7 ) (v x ,y(t) - S x ,y(t)) + W x ,y(t) (26) 

where the friction coefficient 7 is estimated to be 7 ~ 1 day -1 , S x ,y(t) stands for the seasonal 
additive component of the wind velocity and w Xty (t) is the error (noise) term. Whatever the 
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FIG. 4: Semi-logarithmic representation of standardized velocity increment pdf at various scales. 
From top to bottom one goes from small to coarse time scales. All graphs have been vertically 
shifted for sake of clarity, (a) Logarithm of the pdf of the wind velocity increments 5 T v x for the 
Eindhoven wind series. Time scales r go from 1 hour to 5 days, (b) Same plots as in (a) but for the 
increments of the longitudinal velocity field 5 T vu in a high Reynolds number turbulence experiment 
(see text). Scales r extend from the Kolmogorov scale to the integral scale. 

precise definition used to compute the local fluctuations, the results presented in the next 
section remain unchanged. 



B. Evidences for a mesoscale cascade 



As recalled in the introduction, unlike inertial subrange turbulence, few papers have been 
devoted to the study of scaling properties of atmospheric fields from moderate to large scales. 
In Ref. [4J, the authors find that cloud satellite data display multiscaling over the range 
1 — 5.10 3 km and suggest the existence of a (anisotropic) cascade process from planetary 
scales down to the microscales. As far as scaling properties of surface layer wind speed 
are concerned, some studies focus on (multi-) scaling properties. In Ref. [7|, the authors 
use Detrented Fluctuation Analysis on various hourly averaged wind series and provide 
evidences of a cross-over between two scaling regimes separating mesoscale range from very 
large scale range. The scale of the cross-over was found to be T ~ 5 days and the author 



conjectured the possibility of the existence of multiscaling for scales below this scale, 
case the scale T could be identified to some integral (injection) scale. In the Refs. 



'n that 
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Lauren et al. perform an analysis and a modeling of atmospheric wind in both mesoscale and 
microscale (turbulent) regimes. They show that the concepts of multiscaling and cascades 
are also pertinent for characterizing and simulating low-wavenumbers properties of surface 
layer winds. 

Though the above mentionned studies provide hints for the existence of intermittency in 
wind fluctuations at large scales, as explained in section [Til because of diurnal oscillations, 
discreteness of the data, limited size of the scaling range and for purely statistical consid- 
erations, the reported estimations of (multi-) scaling properties of mesoscale wind data are 
far from being reliable. A more direct and convincing illustration of the intermittent nature 
of mesoscale wind variations is provided in Fig. H^a) where we have plotted the normal- 
ized pdf of the increments of (deseasonalized) wind components v x of the Eindhoven site 
in logaritmic scale. One sees that increment distributions go from 'stretched exponential" 
like functions (with large tails) to Gaussian like behavior (parabola) when going from 1 
hour to few days time scales. This behavior is very similar to the features observed in fully 
developed turbulence laboratory experiments. In Fig. IU(b), for comparison purpose, we 
have plotted in the same logarithmic scale, the standardized pdf of longitudinal velocity 
increments computed from experimental data obtained by B. Chabaud and B. Castaing in 
a low temperature gaseous Helium jet experiment 52] (the Taylor scale based Reynolds 
number is R\ = 929). As the scale r varies from the dissipation to the integral scale, one 
observes the same deformation of the pdf from large tailed shape to Gaussian-like shape. 
As explained is section Hi Al (Eq. ( flUj) ). such a variation of the pdf behavior across scales 
is usually associated with the existence of an intermittent cascade. Let us notice that the 
smallest time scales we considered, correspond to scales greater than or equal to the "injec- 
tion" largest scale of atmospheric boundary layer turbulence. Consequently, according to 
our observations, velocity increment statistics at large (atmospheric surface layer) turbulent 
scales are characterized by a strong kurtosis. This contrasts with the situation in laboratory 
experiments where the distribution are nearly gaussian at large scales (as illustrated by the 

n n n 

bottom graph in Fig. IU(b)). Similar observations have been performed in [8j, |9j. In [8(, 
the authors interpret the intermittency of large scale atmospheric turbulence by the fluc- 
tuations of the turbulence intensity at this scale. The velocity mean is no longer constant 
but stochastic. Similar observations on the vorticity field have been made in . In the 
following, we suggest that these fluctuations of turbulence intensity can be interpreted as 



18 



the result of a cascading process starting at a larger time scale. 

As advocated in section Hi C\ the best way to reveal the presence of an underlying random 
cascade and to estimate the intermittency coefficient is to study the magnitude correlations 
functions associated with velocity small scale variations. If one writes equation (|14|) for the 
small scale increments of respectively v x and v y (where the seasonal components have been 
removed) one can define two magnitude processes and two processes Q x and Q y : 

U x , T (t) = H\5 T V x (t)\) = n x , T (t)+ H\e x (t) I) 

u VtT (t) = H\5 T v y (t)\) = n y , T (t) + H\e y (t)\) 

If one assumes that the noises e x and e y are independent gaussian white noises, then thanks 
to the fact that Var(ln(|e|)) ~ 1.23, one can compute the correlation coefficient of Q x and Q y 
from the correlation of u x and uo y . The obtained results for all data series are summarized 
in table [Til One clearly sees that the correlation coefficient of magnitude processes is for all 
series close to 20 % while the estimated coefficients for the Q components are close to 1 2 . 
This result strongly suggests that the processes Q x and Q y are identical and therefore Q is 
a scalar quantity. One thus has (up to season components): 

S T v x (t) = e n ^e x (t) 
5 T v y (t) = e n ^e y (t) . 

The scalar c Qt ^ is simply the (stochastic) amplitude of velocity fluctuation vector at scale 
r. 

According to these findings, in the following we compute a surrogate of the scalar process 
Q(t) as follows: 

Or(t) = X - In (6 T v x (t) 2 + 6 T v y (t) 2 ) (27) 

which is less noisy than individual magnitudes u x or u y . Let us notice, as already mentioned, 
that if one replaces in previous analysis, S T v x ^ y by the error term in a Langevin model of 
y XiV as in Eq. ( |26l) or by small scale wavelet coefficients, all the results remain unchanged 



501 ] . Since the results presented below do not depend on the chosen small scale r (we focus 



2 notice that estimated coefficients greater than 1 are probably due to our hypothesis concerning the nor- 
mality of e X)J/ ; such an assumption can only be a rough approximation for some series because of finite 
size and granularity effects 
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Location 


lo correlations 


f! correlations 


Vignola (Ajaccio) 


0.19 


1.16 


Ajaccio 


0.22 


1.31 


Bastia 


0.20 


1.10 


Calvi 


0.14 


1.10 


Conca 


0.22 


1.07 


Figari 


0.17 


1.32 


Renno 


0.21 


0.71 


Sampolo 


0.11 


0.68 


Eindhoven 


0.22 


1.01 


Ijmuiden 


0.22 


0.99 


Schipol 


0.21 


0.98 



TABLE II: Correlation coefficients of magnitude components for the different sites 



on properties involving mainly lags greater than r), we omit any reference to this scale and 
denote Q(t) the local "magnitude" estimated at small scale according to Eq. (1271) . 

The process Q(t) also possesses a seasonal component which means that seasonal effects 
manifest not only through an additive (locally) periodic component but also through a diur- 
nal modulation of the velocity stochastic amplitude. This multiplicative seasonal component 
has been estimated using the same local least squared method as for the additive compo- 
nents 50J. In the following, we assume that all seasonal effects have been removed from the 
estimated field Q(t). 

Along the line of the method described in section III CI we have estimated the magnitude 
Q(t) correlation functions p(At) = Cov(f2(t), Q(t + At)) for all the time series. As illustrated 
in Fig. [51 when one plots p 1 ^ 2 (At) as a function of ln(At) one observes, for all series, a 
decreasing linear function that becomes zero above some lag T. It is striking to observe that 
all slopes are close to each other and that the "correlation" scale (integral scale) is roughly 
the same for all sites. In order to handle less noisy curves, we have plotted in Fig. [6](b), 
within the same representation, the mean correlation function for all sites in Corsica, for 
all sites in Netherlands and for the high frequency series at Vignola (in Fig. E](a) the same 
graphs are displayed using a linear representation). Up the some remaining bump around 
the lag of one day due to the non perfect removing of the seasonal components, one observes 
a well defined linear dependence on a range [ln(r), ln(T)]. This means that for all wind 
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2 4 

ln(At) (hours) 



FIG. 5: Square root of magnitude covariance functions estimated for all wind series. The three 
bottom graphs correspond to Netherlands hourly series while the top graph is correlation of magni- 
tudes associated with the 'high frequency' Vignola data series. All the graphs have been arbitrary 
shifted vertically for clarity purpose. One can observe comparable values of the parameters /3 2 and 
T (see text). 



series, the covariance of fl(t) reads (At > r): 

p(At) = /3 2 ln 2 



. T , 
•At' 



In section III D\ we have explained how such a square logarithmic dependence of the single 
point covariance can be the result of log-correlated Eulerian and Lagrangian fields. This 
feature has also been observed on various laboratory turbulence data 12|, |40j. These ob- 



servations are therefore direct evidences that a random cascade mechanism can be involved 
in the energy transfer at scales much greater than the usual turbulent microscale. Let us 
remark that the value of the integral scale we found is T ~ 5 days in Holland and Corsica 
and the value of the slope (3 ~ 0.07. In order to deduce Eulerian (or Lagrangian) value of 
the intermittency coefficients, one would need to know the value of — where r\ is a small 
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FIG. 6: Mean magnitude correlation functions associated with wind series in Corsica (o) and 
Netherland (•). The solid line represents the magnitude covariance of the Vignola series (10 min 
rate), (a) The graphs are in linear scales, (b) The square root of the covariances are represented 
as function of the logarithm of the time lag At. 

time scale cut-off above which Lagrangian variations are no longer intermittent. For exam- 
ple if one sets T ~ 5 days and one chooses rj to be a typical time scale separating micro 
and meso scales, i.e., rj ~ 10 min, according to Eq. one gets an Eulerian intermittency 
coefficient A 2 ~ 0.03, i.e., very close to the value found for fully developed turbulence (see 
e.g., Q,Q,Q)- 



IV. CONCLUSION AND PROSPECTS 



The goal of this paper was to provide evidences that the surface layer wind fluctuation 
statistics in the mesoscale range are, very much like microscale turbulent statistics, charac- 
terized by strongly intermittent properties. We have first reviewed how the intermittency 
of continuous random cascades can be advantageously studied using magnitude correlation 
analysis as compared to standard scaling or magnitude cumulant estimations. Within this 
framework and using various wind velocity and direction time series in Corsica (France) and 
Netherlands, we have shown that one can define a scalar magnitude field that displays " uni- 
versal" squared log correlation functions. Such peculiar shape of time correlation functions 
at a fixed position is shown to result from a continuous cascade (with log-correlated mag- 
nitude) that governs the fluctuations in both Eulerian and Lagrangian frames. The same 
behavior has been observed in laboratory turbulence experiments. The existence of some 
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mesoscale "energy cascade" and other similarities with the 3D isotropic turbulence proper- 
ties, raises various fundamental questions. Unlike synoptic circulations, mesoscale motions 
can be associated with a wide variety of phenomena covering a large range of characteristic 
scales from thunderstorms, mountain waves to front dynamics. In that respect, the "univer- 
sal" features discussed in this paper have to be confirmed using further experimental data. 
Let us note however that a characteristic time scale of few days is usually associated with 
front dynamics fl and some comparable typica. time scales have been already reported in 
the literature |5|, l7(. From a practical point of view our findings lead to a better charac- 
terization of the statistical properties of wind fluctuations. The framework of intermittent 
statistics and random cascade models can be applied to address problems related to wind 
resources assessing, extreme events characterization or to design a simple stochastic model 
in order to perform short term wind predictions 



50. 
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